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Abstract 

We have studied second harmonic generation (SHG) in a polar ferrimagnet GaFeC>3, employing a 
FeOg cluster model in which the Fe atom is slightly shifted from the center of the octahedron. The 
electric-dipole transition could take place between the 3d states through the effective hybridization 
of the Ap states with the 3d states, due to the breaking of the space-inversion symmetry In the 
third-order perturbation with H m t = — -j • A, we calculate the probability per unit time, I^aa, f° r 
the process that two photons are absorbed with polarization parallel to the a axis and one photon 
is emitted with polarization parallel to the n (= a, b, c) axis. The calculated SHG intensities consist 
of several peaks as a function of two-photon energy in agreement with the experiments. It is found 
that the corresponding amplitude S aaa at each Fe site changes its sign while Sbaa remains the same 
with the reversal of the direction of the local magnetic moment. This implies that I aaa would 



disappear while Ib aa would survive in the paramagnetic phase in accordance with the experiment. 
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I. INTRODUCTION 

When the space-inversion symmetry and the time-reversal symmetry are simultaneously 
broken, novel magneto-optical effects are expected to come out in the optical spectroscopy. 
Such effects are known as the Kerr effect, the Faraday effect, the reciprocal dichroism, the 
magneto-chiral dichroism, and so on.-~- 

Since GaFeC>3 exhibits simultaneously spontaneous electric polarization and magnetiza- 
tion at low temperatures, this compound is quite suitable to investigate such magneto-electric 
effects. Remeika was the first who synthesized the compound for decades ago. 6 Rado ob- 
served the large magneto-electric effect.- Recently, untwinned large single crystals have been 
prepared. 8 The optical absorption measurement has been carried out with changing the di- 
rection of magnetization,- and the magneto-electric effects on the absorption coefficient have 
been observed. That is, the absorption spectra depend on the direction of the magnetization 
in the region of photon energy 1.0 — 2.5 eV. In our previous paper,— we have analyzed the 
spectra through the microscopic calculation. The calculated spectral shape has agreed well 
with the experimental curve.- 

In the field of the nonlinear optics, it is known that the breaking of the spatial in- 
version symmetry in polar or chiral materials gives rise to second harmonic generation 
(SHG).— ~— When the time-reversal symmetry is simultaneously broken in polar magnetic 
materials, the SHG intensities with some specific polarization are activated by the presence 
of magnetization.—"— In this context, the SHG spectra have been observed and analyzed in 
plenty of systems such as Cr203 and multiferroic materials,—^— and magnetic field induced 
SHG has been investigated in several semiconductors.— - — Recently, the SHG experiments 
have been carried out in GaFeOa.— ~— Ogawa et al. have measured the magnetization- 
induced SHG spectra of GaFeOs in the region of the two-photon energy 2.5 — 4.5 eV.— The 
purpose of this paper is to elucidate the origin of the SHG spectra through the microscopic 
calculation. 

It is known that the crystal of GaFeOs belongs to the space group Pc2\n and has an 
orthorhombic unit cell. 32 Each Fe atom is octahedrally surrounded by O atoms, and is 
slightly displaced from the center of the octahedron along the b axis. There are "Fel " and 
"Fe2" sites for Fe atoms, where the displacement is 0.26A at Fel sites and — 0.1 1A at Fe2 
sites. 8 Thereby the spontaneous electric polarization is generated along the b axis. The local 



magnetic moments at Fel and Fe2 sites are known to align antiferromagnetically along the 
±c axis. Note that the actual compound deviates slightly from a perfect antiferromagnet 
and behaves as a ferrimagnet.— The origin for this deviation is not fully understood, but 
is inferred that the Fe occupation at Fel and Fe2 sites are slightly different from each 
other.- We assume the system as a perfect antiferromagnet in the following analysis. As 
shown later, we could obtain the SHG intensities in the antiferromagnetic phase, since the 
system is breaking the space-inversion symmetry. This situation is different from those of 
EuTe and EuSe, 34 where no SHG intensities exist in the antiferromagnetic phase because Eu 
atoms reside on centrosymmetric sites. We expect that the slight deviation from the perfect 
antiferromagnet would cause merely minor quantitative change in the SHG intensities. 

We introduce a FeOe cluster model in which the Fe atom is slightly displaced from the 
center of the octahedron. 10 We neglect the slight distortion of octahedron. We take account 
of the Coulomb interaction between the 3d states, the spin-orbit interaction on the 3d states, 
and the hybridization of the oxygen 2p states with the Fe 3d and 4p states. The same cluster 
model has been successfully applied to investigating the magneto-electric spectra in the 
optical absorption^ - and the directional dichroic spectra in the K-edge x-ray absorption^ in 
GaFe03. We evaluate the matrix elements of the electric dipole (El) transition between the 
3d and the 4p states using the atomic Hamiltonian from the conventional form H int = — -j-A 
where j is the current operator and A is the vector potential. As shown in Sec. Ill, the matrix 
elements thus evaluated are found larger than those evaluated from another conventional 
form H mt = — P • E where P is the electric dipole operator and E is the electric field. The 
situation may become different without using the atomic Hamiltonian in the band structure 
calculation. The El transition between the 3d states could eventually take place through the 
effective 4p-3d hybridization due to the breaking of the space-inversion symmetry. Thereby 
the effective El transition matrix elements become larger than those of the magnetic dipole 
(Ml) transition in the 3d 5 configuration. In the present cluster model analysis, we regard 
ifint = — -j ■ A more fundamental. 

Bearing the above situation in mind, we avoid to use the form H mt = — P • E in the 
cluster model analysis. Then, using perturbation theory to third-order with H int = — -j ■ A, 
we formulate the SHG intensity from the probability per unit time of the process that the 
incident two photons with the frequency u are absorbed and one photon with the frequency 
2cu is emitted. Although this approach seems different from the conventional analysis using 



the nonlinear susceptibility, the expression for the probability amplitude is quite close to 
that obtained from the nonlinear susceptibility. The difference is that the El transition 
matrix elements in the formula are not evaluated from H mt = — P ■ E. On the basis of this 
formula on the cluster model, we calculate the SHG intensities as a function of two-photon 
energy on various polarization conditions in the so-called phase-matching condition. The Ml 
transition is found to give a minor contribution to the SHG intensity. The calculated SHG 
spectra show multi-peak structure in the 1.0-4.5 eV range in agreement with the available 
experimental data. Another finding is that, when the polarizations of both the incident and 
emitted photons are all parallel to the a axis, the amplitude (complex number) at each Fe site 
change its sign with the reversal of the direction of the local magnetization. This indicates 
that the SHG intensity, which will be denoted as J aoo , would disappear by passing through 
from the antiferromagnetic phase to the paramagnetic phase. Therefore, this spectrum may 
be called as magneto-electric one because it is activated by the magnetization. On the other 
hand, we confirm that, when the polarization of incident photons is parallel to the a axis 
while that of the emitted photon is parallel to the b axis, the corresponding amplitude at 
each Fe site keeps the same value with the reversal of the direction of the local magnetic 
moment. This implies that the SHG intensity, which will be denoted as It, aa , would change 
little by passing through from the antiferromagnetic phase to the paramagnetic phase. These 
characteristics of the polarization dependence as well as the spectral shape agree with the 
SHG spectra observed in the reflection measurement,— 1 ^ although our results is not for the 
reflection spectra. 

This paper is organized as follows. In Sec. II, we introduce a cluster model FeOg. In Sec. 
Ill, we describe the optical transition operators associated with Fe atoms. In Sec. IV, we 
derive the formula of the SHG intensity, and present the calculated spectra in comparison 
with the experiment. The last section is devoted to concluding remarks. 

II. HAMILTONIAN FOR A FeO e CLUSTER 

Considering a FeOg cluster, here we briefly summarize the corresponding model Hamil- 
tonian. The details are found in our previous papers.— 1 ^ It may be expressed as 

H = H 3d + H 2p + H 4p + Hl% 2p + Hi% 2p . (2.1) 



The H 3d represents the energy of Fe 3d electrons, which includes the intra-atomic Coulomb 
interaction expressed in terms of the Slater integrals, the spin-orbit interaction, and the 
energy arising from the exchange interaction via the exchange field from neighboring Fe 
atoms. The energy of the Fe 4p states and that of the oxygen 2p states are represented by 
H Ap and H 2p , respectively. 

The H h r p and H h p ^ p represent the hybridization energies of the O 2p states with the 
Fe 3d and 4p states, respectively: 

Hl% 2p = E t^ P (j)dl aP ^ + H.c, (2.2) 

jrja-m 

K* 2P = E <C P (MW + H.c , (2.3) 
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where d^ and p'^ a stand for creation operators of electron with the 3d orbital 



m = x 



y 2 , 3z 2 — r 2 , yz, zx, and xy) having spin a and of the local 4p orbital (rj = x, y, and z) having 
spin a, respectively. The pj Va is the annihilation operator of electron with the oxygen 
2p orbital at neighboring site j. The sum over j is taken on neighboring O sites. The 
hybridization parameters t^~ 2p (j) and t p ~ p (j) are expressed in terms of the Slater- Koster 
two-center integrals. 

The Fe atom is slightly displaced from the center of the octahedron; the shift 5 is 0.26A at 
Fel sites and — 0.11A at Fe2 sites along the b axis. Therefore, the hybridization parameters 
are modified. We evaluate the modified Slater-Koster two-center integrals for the Fe atom 
by assuming that (prfcr) 2p ,3d, (pd7i) 2p ,3d ex d~ A , and (ppcr) 4 p )2 p, (pp^ip^p oc d~ 2 for d being 
the Fe-0 distance. 3 - With these modified values, we obtain the ligand field Hamiltonian on 
the 3d states, which is deviated from the cubic symmetry. In the second-order perturbation, 
it may be given by 

H 3d ~ 3d = £ &dLd m > ff , (2.4) 
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with 



W = E C 2P 0')^ 2P 0')/A, (2.5) 
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where A denotes the charge transfer energy. In addition to the ligand field Hamiltonian, we 
have the effective hybridization between the 4p and 3d states, due to the breaking of the 
space-inversion symmetry. In the second-order perturbation, the hybridization energy may 



be given by 

^ 4p - M =E^ 3 Xv^ + H.c, (2.6) 

rfrau 

with 

T4p-3d _ ST^ L v'y \J) L mr) \J ) ,_ _s 

rfm ~ 2-^4 £]4p _ ^2p ' \ •') 

m 
where E 4p and E 2p are the average of the 4p-band energy and the energy of the 2p electron, 

respectively. The denominator in Eq. (12. 7p is approximately estimated as E 4p — E 2p « 17 
eV. 

In the following numerical calculation, we use the same parameter values as in our pre- 
vious papers,— 1 ^ except for A = 4.0 eV, which is slightly larger than the previous value 3.3 

eV. 



III. INTERACTION BETWEEN ELECTROMAGNETIC WAVE AND ELECTRON 

We concentrate our attention on Fe atoms. Then, the interaction between electrons and 
the electromagnetic wave with polarization vector e and wave vector q is approximated as 

#int = -- Y] j(q, i) ■ A(q, i) + H.c, (3.1) 

i 

with 

a ] n {i)a n ,{i), (3.2) 



j(q,i) = Y, f^-^nnir-nWr-Ti 



A(q,i) = W^ec q e^ r % (3.3) 

where rj is the position vector of the Fe atom at site i, and a n (i) is the annihilation operator 
of electron with the local wave function 0„(r — r$). The j nn '(r — r,) in Eq. (13.21) is given by 

ten 6 

j nn /(r - Ti) = — -[(V#;)0 n / - 0n V <M - ~ A ^n^' 

+ — cV x [0;S<M, (3-4) 

mc 

where e, m and /iS mean the charge, the mass, and the spin operator of electron, respectively. 
The interaction Hamiltonian can be rewritten as 



tfmt = -e,/-|^^T(q,e,z)c q e jC l^ +H.c, (3.5) 
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where the transition operator T(q, e,i) is defined as he ■ j(q, i)/e. 

First, we consider the El transition. Putting e* q '( r ~ r j) = 1 in Eq. (j3.2p . we evaluate the 
first term in Eq. (13.41) by using the relation 
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,, — w d 3 r = — ^(e„ - e„') / #;z0 n /cl 3 r, (3.6) 



where e n and e n / describe the energy eigenvalues with eigenfunctions <fi n and 4> n i, respectively. 
In the present study, <p n and <f> n > are assigned to the 4p and 3d states. Then, the transition 
operator T E1 is summarized as 

T E \ q ,e,t) = iB El J2 N^\p'l{i)d ma {i) - dLWP^Wl. (3-7) 

with 



/■oo 

5 S1 = (e 4p - e 3(i ) / r 3 J R 4p (r) J R 3 d(r)dr, (3.8) 

Jo 

where R 3c i(r), -R 4p (r) are radial wave-functions of 3d, Ap states with energy e^d, Qp in the 
Fe atom. Within the HF approximation in the ls 2 3d 5 4p 0001 -configuration of an Fe atom,— 
we estimate it as B E1 « 7.7 x 10~ 8 cm-eV. Coefficient NE^s are given by N E ^ 2 _ 2 = l/\/5, 
N E l z 2_ r 2 = — 1/y/TE, A 7 ^ = l/y/E, AT^. = 1/Vo for polarization parallel to the a; axis, 
iV x % = 1/V5, N E y y2 = -l/y/E, N E l z2 ^ r2 = 1/y/TE, N E l z = l/y/E for polarization parallel 
to the y axis, and N E / ZX = l/y/E, N E y Z = l/y/E, N E ^ z2 _ r2 = 2/y/lE for polarization parallel 
to the z axis, respectively. 

The matrix elements of the El transition are sometimes evaluated from the dipole inter- 
action, 

i7int = -^P 4 -E(q,*), (3.9) 

where Pj is the dipole operator for electrons of the Fe atom at site i, and E(q, i) is the electric 
field. In Eq. (13. 9p . the matrix element of the dipole operator between the 4p and 3d states 
may be estimated as Eq. (13. 8p divided by e 4p — e^d, while E(q, i) = — - — ^— = — A(q, i) for 
the oscillating field with the frequency u. Therefore, the use of Eq. (13. 9p underestimates the 
.El-transition matrix elements by a factor fvjj/[e^ v — 63d) with e 4p — 63^ > 10 eV and hu ~ 
several eV's. We think Eq. (13.11) more fundamental from the microscopic standpoint. 

Now we evaluate the matrix elements of the El-transition between the 3d states. The 
matrix elements between the 3d states take finite values when the 3d states mix with the 4p 
states. As a first step of the evaluation, we calculate the energy eigenstates |$ n (<i 5 )) with 



eigenenergy E n (d 5 ) in the 3<i 5 -configuration, and |$ n (<i 4 )) with eigenenergy E n (d 4 ) in the 
3<i 4 -configuration, by diagonalizing the Hamiltonian H^ + H 3d ~ 3d , If the displacement is 
neglected, Fe atoms are under the cubic symmetry. The displacement gives rise to an addi- 
tional trigonal field, which makes the energy levels split further. The spin-orbit interaction 
and the exchange field further modify these states. Note that the matrix elements of the El 
transition do not exist between these states. 

Next, we treat the effective hybridization H 4 P~ 3d within the first order perturbation. The 
modified wave function \^ n (i)) may be written as 

|* n (i)> = |$ n (rf 5 )) 

where |$ m (<i 4 ), kr/cr) represents the state of four electrons in the 3d states and one electron 
in the 4p states specified by rj(— x,y, and z), spin a, and momentum k. The sum over 
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k may be replaced by the integral with the Ap DOS, which is explicitly given in Ref. 
>From Eq. f 1 3 . 1 1) . the matrix elements of the El transition between the states in the 3d 5 - 
configuration are given by 

[T sl (q,e,z)] n> = (* n ,(i)|T B1 (q,e,i)|* n (i)> 

v (^(^ 5 )l^ 1 (q ; e^)|$ m (d 4 ),kr/a)($ m (ci 4 ),k W |g 4 P- M |$ n (d 5 )) 
" ^ E n (d*) - E m (d*) - e 4p (k) 

v (^(^ 5 )l^ 4p - M |^(ci 4 ),kr / a)($ m (rf 4 ),kr /( T|T gl (q,e,z)|$ n (rf 5 )) 
X„ E n ,(dS)-E m (d*)-e 4p (k) 

(3.11) 

Finally, we close this section by evaluating the matrix elements of the Ml transition. In 
doing so, we approximate the third term in Eq. (13.41) as 

/"e iq '< r - ri Vx(^S^)d 3 r = -zqx f foS^e^'^dh 

« -iqx 1 0;S0 n ,d 3 r. (3.12) 

In addition to this term, we have, from the first term of Eq. (I3.4p . the similar term to 
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Eq. (I3.12p . in which S is replaced by L/2 (L is the orbital angular momentum). See Ref. 
for the derivation. Hence, the corresponding transition operator may be expressed as 

T M1 (q,e,i) = i\q\B M1 ^ iV^ mV ,dL« WO, (3-13) 



tram' era' 



where B M1 = h 2 /2m = 3.8 x 10~ 16 cm 2 • eV. When the photon propagates along the 
z axis, N^ m , a , = {ma\ — (L y + 2S y )\m'a') for polarization parallel to the x axis, and 
Nmam'a' = ( m<T \Lx + 2S x \m'a') for polarization parallel to the y axis. The matrix elements 
between the 3d states in the ^-configuration are given by 

[T M1 (q,e,i)] <n = (* n ,(i)\T M1 (ci,e,i)\* n (i)) 

« (<l> n ,(d 5 )\T M1 ( q ,e,t)\<5> n (d 5 )). (3.14) 

These values are found much smaller than the matrix elements of the El transition given 
by Eq. (13~TT|) . 

IV. SECOND HARMONIC GENERATION 

We consider the process that two photons with the frequency oo, wave vector q, and 
polarization e are absorbed and one photon with Co, q, and e is emitted with Co = 2uo, as 
illustrated in Fig. [TJ In the third-order perturbation with H mt = — -j • A, the probability 
per unit time for the process may be expressed as 

2 



I(e,co, q; e, Co, q) oc 



X;^(e,u;,q;e,a;,q;z)e^-^ 



Cd 2 5(Co-2u), (4.1) 



where the amplitude S is given by 



oc 



+ 



S(e, u, q; e, Co, q; i) 

1 s-^f [T*(q,e,i)] g , n >[r(q,e,i)] n / iW [T(q,e,i)] nifl 
wv^T^l ( e «' - e g - Zftu - ^ r )(e n -e g -fwo-iT) 

[T(q, e, i)] gjn , [T*(q, e, i)] n >, n [T(q, e, i)]», 3 

(fn' — e g + frw — fiw — iT)(e n — e g — fkj — iT) 

[T(q, e, z)] 9) „/[T(q, e, i)] n> [T*(q, e, z)] n , s 



(e n / — e g + hu — ftio — iT)(e n — e g + fwo — iT) J 

The T represents the life-time broadening width by other random perturbation on the ma- 
terial system. The e g is the energy of the ground state in the 3d 5 -configuration, and e n is 
the abbreviation of E n (d 5 ). Note that Eq. (14. 2 j) resembles the conventional expression of 
nonlinear susceptibility, which is based on the interaction H[ nt = — ^ P, ■ E(q, i). The sum 
over % in Eq. (14.11) is made on all Fe sites, and the so-called phase-matching condition q = 2q 
has to be satisfied in order to get finite intensities. Note also that, if the wave interaction 



length £ is finite, the momentum conservation would be relaxed within Aq ~ I. In reflection, 
the surface layer with 1/q thickness could contribute significantly to generating reflection 
wave without phase matching-condition.— Since the matrix elements of the El transition 
is much larger than those of the Ml transition as discussed in Sec. Ill, all the transition 
matrix-elements in Eq. (14 .2p could be replaced by those of the El transition. 
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FIG. 1: Schematic description of second harmonic generation. Two photons with e, u), q are 
absorbed, and one photon with e, a), q is emitted. 

We consider the situation that the polarization of incident photons is parallel to the 
a axis. When the polarization of the emitted photon is parallel to the c axis, the SHG 
intensity, which will be denoted as I caa (2u), is found to vanish. On the other hand, when 
the polarization of the emitted photon is parallel to the a axis, non-zero SHG intensity 
-^aaa(2o;) is obtained as shown in Fig. [2J The broken line represents the result with replacing 
T*(q, e, i) by T M1 *(q, e, i) in Eq. (14. 2p . demonstrating that the contribution of the Ml 
transition is much smaller than that of the El transition as estimated in Sec. IIHI The 
spectral shape is composed of multi-peak structure: a small peak around 2u = 1.75 eV, 
a small peak around 2.5 eV, and a two-peak structure around hu = 3 — 4 eV. The two- 
peak structure reasonably captures a whole aspect of the experimental SHG intensity shown 
in the inset of Fig. [21 although the experimental spectrum is obtained on the " S[ n S out " 
configuration in the reflection geometry 28 

For the polarization of the emitted photon parallel to the b axis, we obtain the spectra 
Ibaa (2a;) as shown in Fig. |3j In comparison with the spectral shape of I aaa (2u), the peak 
around 2u = 1.75 eV becomes larger, but the other peaks remain similar to those of I aaa {2uj). 
As shown in the inset of Fig. [31 a large peak is observed around 2a; = 1.5 eV in the SHG 
experiment, 29 which is reproduced by the present theory. Our calculation also implies that 
another peak in the 3.3 — 4.0 eV range is anticipated if corresponding experiment will be 
available. Note that the SHG intensities have been observed around 2.5 — 4 eV on the 

10 



"Sin-Pout" configuration in the reflection experiment.— ^i The obtained intensity is about 
1/4 of I aaa (2u) in this energy region, corresponding well to the ratio of intensity on the 
"Sin-Pout" configuration to that on the " Si n S out " configuration in the reflection experiment. 
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FIG. 2: The second harmonic generation (SHG) intensity I aaa (2u)) as a function of two-photon 
energy 2oj with e and e parallel to the a axis. Inset: experimental SHG intensity in the reflection 
geometry; the s-polarized light is radiated with incident ang le 5 degree on the surface of ac plane, 



and the s-polarized reflected SHG light is observed (Rcf. 



M 



Now we examine what happens to the SHG intensity when the system undergoes a phase 
transition into the paramagnetic phase. The probability amplitude S aaa (2u]i) is found to 
change its sign, when the direction of the local magnetic moment is reversed. This indicates 
that 'Y^ li S aaa {2oj\i) would be canceled out when the local magnetic moment is randomly 
oriented, and that I aaa (2u) would disappear in the paramagnetic phase, in agreement with 
the experiment. For this reason, I aaa (2u) may be called as the magnetization-induced SHG 
intensity, which corresponds well to the spectra around 2hu ~ 3 — 4 eV for the " SinSout " 
configuration in the reflection experiments—^- On the other hand, the amplitude Sb aa (2u; i) 
is found to remain the same when the direction of the local magnetic moment is reversed. 
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FIG. 3: The SHG intensity I^ aa (2oj) as a function of two-photon energy 2a; with e and e parallel 
to the a and b axes, respectively. Inset: experimental SHG intensity (Ref. |29|). 



This leads to that Y2i <Sftaa (2cu; i) would not be canceled out when the local magnetic moment 
is randomly oriented, and that Ib aa (2oj) remains finite in the paramagnetic phase, which 
corresponds well to the spectra around 2hw ~ 3 — 4 eV for the ".Sin-Pout" configuration in the 
reflection experiments.— 1 ^ 

Finally we comment on the energy levels. In the present calculation scheme, if we disre- 
gard the displacement of Fe atoms and the spin-orbit interaction, we have the ground state 
characterized as 6 A 1; and excited states characterized as 4 T 1; 2 T 2 and 4 T 2 with excitation 
energies 2.09, 2.16 and 2.8 eV, respectively. These states, however, have no El transition 
matrix elements from the ground state. Only after taking account of the displacement of Fe 
atoms, they have finite El transition matrix elements, but the energy levels may be shifted 
and split. To demonstrate this point, we calculate the absorption coefficient I a b s (w) from 
the formula 



J a6s (w) oc -y]y2\[T E1 (q,e,i)] n J 2 5(e n -e g -hu). 



(4.3) 
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Assuming the photon propagates along the c axis, the polarization is summed over e parallel 
to the a and b axes. Figure S] shows the calculated spectrum (solid line). It is clearly seen 
that the peaks are considerably separate from the positions of 4 Ti, 2 T2 and 4 X2 shown by 
the vertical solid lines. Therefore, it does not seem meaningful to assign directly these peaks 
to the fictitious levels which have no El transition matrix elements. These peak positions 
depend on parameters such as Slater integrals and the hybridization between Fe and O 
atoms. The broken line represents the spectrum calculated from a different parameter set 
that F 2 and F 4 are reduced by multiplying a factor 0.81 instead of 0.9 and A is set to be 4.4 
eV instead of 4.0 eV. Vertical broken lines indicate the corresponding energies of fictitious 
levels. In the experiment, 28 a peak is found around 1.6 eV, and a shoulder structure around 
2.0 eV. The calculated second peak may correspond to the shoulder in the experiment. 
Further adjustment of parameter sets may improve the calculated spectra, but we would 
not seek optimal parameter sets in this paper. 




one photon energy [eV] 

FIG. 4: Absorption coefficient as a function of one photon energy. The solid line represents the 
calculated spectrum. Vertical solid lines indicate the energies of the fictitious levels with disre- 
garding the displacement of Fe atoms. The broken line represents the calculated spectrum with a 
different parameter set (see the text). The broken vertical lines indicate the energies of fictitious 
levels corresponding to the latter parameter set. 
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As regards the fiction levels with higher energies, we have A A\ and 4 E both with energy 
3.68 eV, 2 A 2 with 3.70 eV, 2 7\ with 3.75 eV, and so on. As the same as the low energy 
levels mentioned above, these levels have no El transition matrix elements from the ground 
state, and would be shifted and split due to the displacement of Fe atoms and the spin-orbit 
interaction. Therefore, it would be difficult to assign directly these levels to the peaks on 
the SHG spectra. 

V. CONCLUDING REMARKS 

We have analyzed the SHG spectra in a polar ferrimagnet GaFeOs, using the FeOg cluster 
model where the Fe atom is displaced from the center of the octahedron. We have fully taken 
account of the Coulomb interaction between the 3d states, the spin-orbit interaction on the 
3d states, and the hybridization of the oxygen 2p states with the 3<i and 4p states. The 
El matrix elements between the 3d and 4p states are evaluated on the Fe atom with the 
interaction H mXj = — -j • A. They are larger than those evaluated with the conventional 
form H [nt = — P • E. The El matrix elements between the 3d states could become finite 
through the effective hybridization between the 4p and 3d states owing to the breaking of the 
space-inversion symmetry. Note that the same cluster model has been successfully applied 
to analyzing not only the optical absorption^ but also the A-edge x-ray absorption.— 

In the third-order perturbation with H lTLt = — -j • A, we have derived the formula of the 
probability per unit time for the process that two photons are absorbed and one photon is 
emitted. On the basis of this formula, we have calculated the spectra as a function of the 
two-photon energy in the phase-matching condition. The calculated SHG intensities exhibit 
multi-peak structure, which is in accordance with the experiments.—^ 1 ^ The intensities 
also signify that another peak structures will be found outside of the published experimental 
surveys^ 1 ^ 1 ^; one centered around 1.5 — 2.0 eV in I aaa (2u; i) and the other centered around 
3.0 — 4.0 eV in Ib aa (2u,i). We have found that the amplitude S aaa (2u;i) changes its sign 
while St, aa (2u; i) retains the same value, when the direction of the local magnetic moment is 
reversed. This indicates that I aaa {2u}) would vanish but Ii )aa {2u) would remain finite in the 
paramagnetic phase. Hence our results have reproduced the experimental observations and 
I aaa (2u) could be called as the magnetization- induced SHG intensity, which is completely 
governed by the El process since the Ml contribution is negligible. 
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The quantum-mechanical treatment in the present paper has not directly been applied 
to the reflection and the refraction problem, since the phase coherence is not appropriately 
treated. More elaborate treatments using the coherent state might be required. 11 Researches 
along this line are relegated to the future study. On the other hand, the semi-classical 
treatment is known to cope well with the reflection and the refraction through the nonlinear 
susceptibility.— 1 ^ We have not directly used the conventional semi-classical method, since 
the form H int = — P ■ E underestimates considerably the El transition. Nonetheless the 
probability amplitude is very close to the conventional form of the nonlinear susceptibility. 
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